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The conductance of a two-dimensional electron gas at the transition from one quantum Hall 
plateau to the next has sample-specific fluctuations as a function of magnetic field and Fermi energy. 
Here we identify a universal feature of these mesoscopic fluctuations in a Corbino geometry: The 
amplitude of the magnetoconductance oscillations has an /h resonance in the transition region, 
signaling a change in the topological quantum number of the insulating bulk. This resonance 
provides a signed scaling variable for the critical exponent of the phase transition (distinct from 
existing positive definite scaling variables). 
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I. INTRODUCTION 

A two-dimensional electron gas in a strong perpendicu- 
lar magnetic field has an insulating bulk and a conducting 
edge. While the conductance Gedge for transport along 
the edge is quantized in units of e^//i, the conductance 
Gbuik for transport between opposite edges is strongly 
suppressed. This is the regime of the (integer) quantum 
Hall effect. -"^ Upon variation of magnetic field or Fermi en- 
ergy, a transition occurs in which Gedge varies from one 
quantized plateau to the next, accompanied by a peak in 
Gbuik- 

The quantum Hall plateau transition is not smooth, 
but exhibits fluctuations reminiscent of the universal con- 
ductance fluctuations in metals. In small samples the 
fluctuations take the form of sharp peaks, due to resonant 
scattering between opposite edges mediated by quasi- 
bound states in the bulk.^ In larger samples the inter- 
mediate states form a percolating network at the plateau 
transition, resulting in a smooth peak in Gbuik with rapid 
fluctuations superimposed. 

The conductance fluctuations are typically analyzed in 
a context that emphasises their random, sample-specific 
nature.'' Here we wish to point out one feature of these 
fluctuations that has a sample-independent topological 
origin. 

The quantum Hall plateau transition is a topological 
phase transition, because a topological quantum number 
Q changes from one plateau to the next.^^ We show that 
a change in Q is associated with a resonance in the am- 
plitude AG of the magnetoconductance oscillations in a 
ring (Corbino) geometry. Each unit increment of Q cor- 
responds to a resonant amplitude AG = e^//i (times spin 
degeneracies). Our analysis relies on a scattering formula 
for the topological quantum number, which relates Q to 
the winding number of the determinant of the reflection 
matrix. ^^'^^ 

The analytical considerations in Sec. II are supported 
by numerical calculations in Sec. HI. These show, in 
particular, that the plateau transition in a disordered 
system is reentrant: there are multiple increments of Q, 




FIG. 1: Corbino geometry consisting of a ring-shaped two- 
dimensional electron gas (2DEG, grey) in a perpendicular 
magnetic field B. A pair of electrodes (black), at a voltage 
difference V ^ is attached to the inner and outer edges of the 
ring. 

alternatingly -hi and —1, each associated with an /h 
resonance in AG. We conclude in Sec. IV by showing 
that the critical exponent of the phase transition can be 
extracted from the sample-size dependence of the width 
of the resonance, in a way which preserves information 
on which side of the phase transition one is located. The 
Appendices contain details of the numerical calculations. 



II. TOPOLOGICAL QUANTUM NUMBER AND 
CONDUCTANCE RESONANCE 

To isolate Gbuik from Gedge we consider a Corbino 
geometry, see Fig. 1, consisting of a ring-shaped two- 
dimensional electron gas (2DEG) connected to an elec- 
trode at the inner and outer perimeters. The two- 
terminal conductance G = Gbuik is then fully due to the 
current / through the bulk, without contributions from 
the currents circulating along the edges. 

We assume that the outer electrode, at a voltage V 
relative to the inner electrode, is connected via an A^"- 
mode lead to the 2DEG. (For simplicity we do not include 
spin as a separate degree of freedom.) The reflection 
amplitudes Vmn from mode n to mode m are contained in 
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the NxN reflection matrix r. The conductance G = I /V 
(in units of the conductance quantum Go = e^/h) fohows 
from 

G = GoTr(l-rrt) = Go^(l-i?„), (1) 

n 

with U the NxN unit matrix. The reflection eigenval- 
ues Rn G [0, 1] are eigenvalues of the Hermit ian matrix 
product rr^. Away from the plateau transition all RnS 
are close to unity and G <C e^//i. 

The topological quantum number Q G Z of the quan- 
tum Hall effect in a translationally invariant system is 
the Chern number of the bands below the Fermi level. 
An alternative formulation, ^^'^^ applicable also to a flnite 
disordered system, obtains Q as the winding number of 
the determinant of the reflection matrix, 

where ^ = (j) x h/e is the magnetic flux enclosed by the 
inner perimeter of the ring. To deflne the winding num- 
ber one needs to vary (j) at constant magnetic fleld in the 
2DEG, so that r{(j)) is 27r-periodic. The conductance (1) 
then oscillates periodically as a function of ^, with ampli- 
tude AG. (We will examine in the next Section to what 
extent this applies to the realistic situation of a uniform 
magnetic fleld.) 

If the ^-dependence of r is continued analytically to 
arbitrary complex z = e*"^, the winding number 

Q — '^zero ~ '^pole (3) 

equals the difference of the number of zeros and poles of 
detr inside the unit circle \z\ = 1. A pole may annihilate 
a zero, but the difference Uzero — ^poie can only change 
when a pole or zero crosses the unit circle. Unitarity of 
the scattering matrix requires |detr(^)| < 1 for real ^, 
so poles cannot cross the unit circle. A topological phase 
transition, corresponding to a change in Q, is therefore 
signaled by a zero of detr crossing the unit circle, see 
Fig. 2. 

We conclude that at the transition point detr(^) = 
for some real value of ^, which implies that (at least) 
one of the RnS vanishes. If the zeros are well separated, 
and all other R^s remain close to unity, the conductance 
oscillation amplitude AG would thus show a peak oie^ /h 
whenever Q changes by ±1. This transport signature of 
a topological phase transition is the analogue for a 2D 
system of the ID signature of Ref. 14. 

III. NUMERICAL SIMULATIONS IN A 
DISORDERED SYSTEM 

We have searched for the predicted conductance res- 
onances in a numerical simulation. The Corbino disc is 
discretized on a square lattice (lattice constant a, nearest- 
neigbor hopping energy t). Electrostatic disorder is in- 
troduced as an uncorrelated random on-site potential. 
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FIG. 2: Example of a path taken in the complex plane z — e'^^ 
by a zero (blue dot) of detr, as the system is driven through 
a phase transition in which Q ^ Q — 1- At the transition, the 
zero crosses the unit circle and a reflection eigenvalue vanishes 
for some flux (j)o . The speed at which the zero crosses the unit 
circle determines the critical exponent of the quantum Hall 
phase transition, see Sec. IV. 



drawn uniformly from the interval (— eo,eo). A uni- 
form perpendicular magnetic fleld B is introduced by the 
Peierls substitution, with / = Bo?e/h the flux through 
a unit cell in units of the flux quantum. 

Referring to Fig. 1 we take parameters Rout — 200 a, 
i?in = 150a, W = 60 a. For / = 0.005 we have a magnetic 
length Im = {h/eSY^'^ = 5.6 a, sufficiently small that 
the edge states at opposite edges do not overlap. These 
parameter values are representative for a /im-size GaAs 
Corbino disc in a fleld of 5 T. The disorder strength was 
flxed at eo = 0.2 1. 

We calculate the reflection matrix r at energy /i as a 
function of the flux / using the recursive Green func- 
tion technique. This gives the conductance G via Eq. 
(1). The topological quantum number Q is calculated 
from Eq. (3), by locating the zeros and poles of detr 
in the complex flux plane using the method of Ref. 13. 
(We summarize this method in App. A.) Results are 
shown in Fig. 3. Since we are interested in mesoscopic, 
sample-speciflc effects, this is data for a single disorder 
realization, without ensemble averaging. 

The left panels show that we can access regions with 
different topological quantum number by varying flux 
and Fermi energy. This represents a wide sweep in pa- 
rameter space, with plateau regions of constant Q = 
1, 2, 3, 4, 5 separated by transition regions of fluctuating 
Q. The conductance G, which in the Corbino geometry 
is the bulk conductance, vanishes on the plateau regions 
and fluctuates in the transition regions. 

The right panels zoom in on the Q = to 1 transi- 
tion. The conductance oscillations as a function of the 
flux are only approximately periodic, because the mag- 
netic fleld changes uniformly (and not just inside the in- 
ner perimeter of the ring). Still, a dominant periodicity 
A/ = 1.1 • 10^ can be extracted, corresponding to an h/e 
flux through a disc of radius i?eff ^ 170 a, about halfway 






FIG. 3: Topological quantum number (top panels) and conductance (bottom panels) of a disordered 2DEG in the Corbino 
disk geometry. The left panels give a broad range of flux / and Fermi energy /i, while the right panels show a close-up of the 
2 = 0^1 transition. (Note the different color scales for the left and right panels.) The horizontal dashed line indicates the 
three reentrant phase transitions (see Fig. 7). 



between the inner and outer radii of the ring. (For a /im- 
size Corbino disc this would correspond to a field scale 
of AB ^ 10 mT.) 

To search for the predicted signature of the topologi- 
cal phase transition, we calculate the amplitude AG = 
Gmax — Gmin of the couductaucc oscillations over one pe- 
riod (/, / + A/). This amplitude is plotted in Fig. 4 as 
a function of Fermi energy, for different fluxes /. The 
2 = 0^1 transition should give at least one e'^/h peak 
value of AG, which is indeed what is observed. For the 
larger field values there is only a single e'^/h peak, while 
for the lower field values three peaks develop, consistent 
with the 2 = 0^1^0^1 reentrant transition seen 
in Fig. 3 (upper right panel, dashed line). 



IV. DISCUSSION AND RELATION TO THE 
CRITICAL EXPONENT 

We have shown that the sample-specific conductance 
fluctuations at the quantum Hall plateau transition con- 
tain a universal feature of a topological origin: The mag- 



netoconductance oscillations in a Corbino disc geometry 
have an e'^ /h resonant amplitude whenever the topologi- 
cal quantum number Q is incremented by ±1. The strict 
h/e flux periodicity of the magnetoconductance is bro- 
ken in a realistic setting by the penetration of the flux 
through the conducting region, but we have shown by nu- 
merical simulations that the resonances are still clearly 
observable. 

The theory presented here may motivate an experi- 
mental search for the reentrant phase transition that we 
observed in our simulations. Existing experiments^ in a 
GaAs Corbino disc show a maximal amplitude of about 
0.1 /h in the Q = 1 ^ 2 plateau transition, an order of 
magnitude below the predicted value. This may partly 
be due to limited phase coherence, and partly to the self- 
averaging effect of overlapping resonances. In the exper- 
imental geometry the entire inner and outer perimeters 
of the Corbino disc are contacted to current source and 
drain, while our geometry (Fig. 1) has two narrow con- 
tacts — which helps to isolate the resonances and make 
them more easily observable. 

From a theoretical perspective, the evolution of zeros 
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FIG. 4: Dependence of the amplitude of the magnetocon- 
ductance oscillations on the Fermi energy for different mag- 
netic fields. The n-th trace from the bottom shows the 
/i-dependence of AG = Cmax — Cmin in the flux interval 
(/n,/n + A/), with = 490 • 10-' + 2(n - 1)A/ and 
A/ = 1.1 • 10^. Successive traces are displaced vertically, with 
the Gq = /h peak value indicated by a horizontal black line. 

of the reflection matrix determinant, see Figs. 2 and 7, 
provides a new way to analyse the plateau transition in a 
disordered mesoscopic system. In particular, the critical 
exponent v of the quantum Hall phase transition can be 
extracted from the speed at which a zero of detr(2;) 
crosses the unit circle in the complex flux plane, 

lim ^ In |zo| = const x(/i-/ic)W'^/^ + (9(/i-/ic)^. (4) 

Here /ic is the value of the Fermi energy ji at the phase 
transition, where zq = e^^ for some real flux (j). The 
Corbino disc has inner perimeter of radius L and outer 
perimeter of radius L + W ^ with the limit L ^ oo taken 
at constant W. 

We note the following conceptual difference with the 
usual MacKinnon-Kramer scaling: ^^'^^ The MacKinnon- 
Kramer scaling variable is a Lyapunov exponent, which 
is a nonnegative quantity. Our scaling variable In 1 2:0 1 
changes sign at the phase transition, so it contains infor- 



mation on which side of the transition one is located. 

In App. B we present a calculation of the critical ex- 
ponent along these lines. It remains to be seen whether 
this alternative numerical method has advantages over 
the conventional method based on MacKinnon-Kramer 
scaling. It is also still an open question whether the evo- 
lution of the zeros of det r is more accessible to analytical 
methods than the evolution of Lyapunov exponents. 
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Appendix A: Calculation of the topological quantum 
number 

We summarize the method used to calculate the topo- 
logical quantum number Q of the Gorbino disc, following 
Ref. 13. We insert a flux tube $ = (j)h/e at the center, 
without changing the magnetic field in the 2DEG, and 
obtain Q as the winding number of detr. The integra- 
tion (2) over (j) can be avoided by going to the complex 
z = e^^ plane and counting the number of zeros and poles 
of det r in the unit circle. The difference nzero — ^poie then 
directly gives Q. 



1. Analytic continuation to complex flux 

To perform the analytic continuation to complex flux 
we cut the disc as indicated in Fig. 5. We attach virtual 
leads (see Fig. 6) at the two sides of the cut and construct 
the four-terminal unitary scattering matrix 

g ^ ^/37 | ^ ^p^-y^^ 

t'ya ^7/3 "^77 I 
\t5a ts(3 ts^ rssJ 

We choose a gauge such that the entire phase increment 
(j) from the flux tube is accumulated at the cut, between 
terminals 7 and 3. The reflection matrix r of the original 
two-terminal system (from terminal a back to a) is then 
obtained from 

r = rc,a- (A2) 

/, f A f ^77 ^7<5 ~ ^\ ft-fa\ 
[ta^ taS) y^^^ _ J y^^J ■ 

This expression is in a form suitable for analytic contin- 
uation to arbitrary complex z = e^^ . 



FIG. 5: Corbino disc in a uniform magnetic field con- 
taining additionally a flux tube $ inside the inner perimeter. 
We choose a gauge such that the phase increment = e^/h 
due to the flux tube is accumulated upon crossing a cut in- 
dicated by the dashed line. Virtual leads 7, 6 are attached 
at the two sides of the cut. Together with the two physi- 
cal leads a, /3 they define the four-terminal scattering matrix 
(Al). The reflection matrix r of the original two-terminal 
system is obtained from Eq. (A2), in a formulation which can 
be analytically continued to arbitrary complex z = e^^ . 
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FIG. 7: Positions of the zeros (red) and poles (green, clus- 
tered near the origin) of detr, along the path indicated by 
the dashed line in Fig. 3. Arrows point in the direction of in- 
creasing Fermi energy. The numbers 1,2,3 indicate the three 
phase transitions in which the topological quantum number 
Q changes first from to 1, then from 1 back to 0, and once 
again from to 1. 



This results in the expression 




FIG. 6: Construction of a virtual lead, connecting two sides 
of a tight-binding lattice. The virtual lead consists of paral- 
lel ID chains with on-site energy /x and hopping amplitudes 
alternating between +1 and —1. The tight-binding equations 
enforce that every other site in the chain (colored red or blue) 
has the same wave amplitude at energy /i. The wave ampli- 
tudes ipn, V^n at the left and right end of the n-th chain are 
therefore unaffected by the insertion of the virtual lead. We 
use this device in Fig. 5 to convert a two-terminal geome- 
try into a four-terminal geometry, by cutting the lead at the 
center and connecting the ends to terminals 7, 6. 



det r = det ts^ 




det(, 

tsj - l/z rss 



(A4) 



The zeros and poles of det r are therefore those values 
of z for which, respectively, the determinant in the nu- 
merator or denominator of Eq. (A4) vanishes. By equat- 
ing each of these determinants to zero we obtain two 
generalized eigenvalue problems, with the zeros given by 
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and the poles by 
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2. Evolution of zeros and poles of detr 



The determinant of the reflection matrix (A2) is com- 
puted as a function of the complex variable z by means 
of the relation 



The numerical solution of generalized eigenvalue prob- 
lems is quick and accurate. As an example, we show in 
Fig. 7 the motion of the zeros and poles along the reen- 
trant phase transition of Fig. 3. 



Appendix B: Calculation of the critical exponent 



det{D - CA-^B) = det (^^ [det A]-^ (A3) 



In Sec. IV we outlined how the critical exponent of the 
quantum Hall phase transition can be extracted from the 
speed at which a zero zq of det r{z) crosses the unit circle 



6 



w 




FIG. 8: Geometry of the network model used to calculate 
the critical exponent. This rolled up strip is topologically 
equivalent to the Corbino geometry used in the main text. 
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FIG. 9: Scaling variable as obtained numerically (data points) 
and as fitted to Eq. (Bl) (black lines). 



in the complex z = e^^ plane. Here we demonstrate that 
this approach produces results consistent with earlier cal- 
culations based on the scaling of Lyapunov exponents. 

To make contact with those earlier calculations, we use 
the same Chalker-Coddington network model^^'^^ (rather 
than the tight-binding model used in the main text). The 
parameter that controls the plateau transition in the net- 



work model is the mixing angle a of the scattering phase 
shifts at the nodes of the network. The transition is at 
Qfc = 7r/4. Disorder is introduced by means of random, 
uncorrelated phase shifts on the links between nodes, 
sampled uniformly from [0,27r). 

The network has a strip geometry, of longitudinal di- 
mension L and transverse dimension W (see Fig. 8). The 
longitudinal dimension has twisted boundary conditions, 
ip e^^ip upon translation over L. The transverse di- 
mension is connected to a source and drain reservoir. 
This geometry is equivalent to a Corbino disc, enclos- 
ing a flux ^ = (j)h/e and with source and drain contacts 
extending along the entire inner and outer perimeter. 

Taking the self- averaging limit L ^ oo of Eq. (4) is 
impractical, but sufficient convergence is reached for an 
aspect ratio L/W = 5 and an average of ln|zo| over 
1000 disorder realizations. Results for the W and a- 
dependence of the scaling parameter A = (W/L) {In \ zo\) 
are shown in Fig. 9. 

We extract the critical exponent u using the procedure 
of Ref. 16, by fitting the data to the scaling law 



p=i 



n 

E 

q=l 



(Bl) 



We have checked that corrections to scaling are insignif- 
icant in our parameter range. The symmetry of the net- 
work model (a time-reversal operation followed by trans- 
lation over half a unit cell) implies that only odd powers 
of a — ac appear in the series expansion {cp = = 
for p even). We truncated the expansions such as to 
keep terms of order {a — ac)^ and lower. The result 
u = 2.56 ± 0.03 is consistent with the value u = 2.59 
obtained in Ref. 16. 
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